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Abstract. We study the separability of the Neumann-Rosochatius 
system on the n-dimensional sphere using the geometry of bi- 
Hamiltonian manifolds. Its well-known separation variables are 
recovered by means of a separability condition relating the Hamil- 
tonian with a suitable (1, 1) tensor field on the sphere. This also 
allows us to iteratively construct the integrals of motion of the 
system. 



The Neumann system is among the widest known and best studied 
integrable systems in Mathematical Physics. It decribes the dynamics 
of a point particle constrained to move on the sphere S"", under the 
influence of a quadratic potential V{x) = | 7^ aj. In 

1859, Carl Neumann [201 showed that the equations of motion of the 
"physical" n = 2 case could be solved using the Jacobi theory of sep- 
aration of variables. It was noticed by Rosochatius (see j22|) that a 
potential given by the sum (with nonnegative weights) of the inverses 
of the squares of the (Cartesian) coordinates can be added without los- 
ing the separability property. The system so obtained is customarily 
called the Neumann-Rosochatius (NR) system. 

More than one century later this separability result was generalized 
to the arbitrary n case by Moser ^H]- The starting point to solve the 
problem was the ingenious introduction of a special set of coordinates 
on 5", called spheroconical (or elliptical spherical) coordinates (already 
used, for = 2, by Neumann). They are defined as follows: For given 
sets of real numbers ai < ai < ■ ■ ■ < a„+i and nonzero xi, . . . , Xn+i, 
the coordinates Xa{x),a = 1, . . . ,n, are the solutions of the equation 
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Later on, it was shown that the NR system could be framed within 
the formahsm of Lax pairs and r- matrices (see, e.g., QIE]). Actually, 
it turns out that introducing the Lax matrix, as a function of the 
Cartesian coordinates Xi,yi, z = 1, . . . , n + 1, by 

(1-2) ^(^) - ( -"^r MA)tV) 

where 
(1-3) 



2 X- ai 2^ \ - ai' 

1=1 1=1 



2^ \ - ai 2 \ ^ \- at 

1=1 \ 1=1 

and identifying the cotangent bundle to 5*" with the submanifold of 
]^2(n+i) (jggj2g(j ^fyy constraints 

n+l n+1 

(1-4) = 1, E^.2/i = , 

1=1 i=l 

the Hamilton equations of motion of the NR system acquire the form 

(1.5) 

where 

As a consequence, the spectral invariants of A^(A) are constants of 
motion. In particular, the quantities 

iri = res| det(iV(A)) , 2 = l,...,n + l 

(known as Uhlenbeck integrals), provide n mutually commuting inte- 
grals of motion that ensure Liouville-Arnol'd integrability of the NR 
system, the physical Hamiltonian being given by 

n+l ^ 71+1 

Hnr = 2 E "^^i + 2 5Z f^'f^r 

i=l i>j=l 

Separation of variables is recovered in this formalism noticing that, on 
S"', the zeroes {Xa}a=i,...,n of the matrix entry /(A) define the sphero- 
conical coordinates, and their conjugate momenta are given (as it will 
be explicitly recalled in Section 0]) by the values of the rational function 
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h{X) for A = Aq. Clearly, each pair of canonical coordinates (Aa, /ia) 
satisfy the separated equation 



In this paper we want to provide a further geometrical interpretation 
of the NR system, based on the notions of bi-Hamiltonian geometry, 
generalizing and refining the approach described in j2]j . We will follow 
a recently introduced set up for the theory of separation of variables 
for the Hamilton- Jacobi equations. In a nutshell, such a framework 
can be described as follows. One considers a symplectic manifold M 
endowed with a (1, 1) tensor field N with vanishing Nijenhuis torsion 
(which we will call an uN manifold, provided that a compatibility con- 
dition between and the symplectic form is satisfied); under suitable 
hypotheses, selects a special subclass of canonical coordinates on M 
(called Darboux-Nijenhuis coordinates) that have the property of diag- 
onalizing N. The condition for the separability of the Hamilton- Jacobi 
equation associated with a Hamiltonian H can be given, according 
to the bi-Hamiltonian theory of separation of variables, the following 
intrinsic formulation. One considers the distribution generated 
iteratively by the action of N on the Hamiltonian vector field Xh as- 
sociated with H, and the two-form d{N*dH). Then Xh is separable 
in the Darboux-Nijenhuis coordinates associated with if and only if 



This scheme, in its basic features, has already been considered in the 
literature HZl HUl CHI HI] and applied to various systems (see, e.g, 
13 El CHI); it is fair to say that, in these papers, the uN mani- 
fold structure is fixed "a priori", and that equation ()1.7|) is seen as a 
condition that selects those Hamiltonians which are separable in the 
"preassigned" Darboux-Nijenhuis coordinates. 

In the present paper we will take a different logical standpoint: We 
will consider a given Hamiltonian H (namely, the NR Hamiltonian) 
and look at (jl.7j) as an equation to determine N (and hence the sep- 
aration coordinates). We shall see that it is indeed possible (and, ac- 
tually, easy) to solve such an equation by means of a couple of natural 
Ansdtze, thus arriving to induce from H the separation coordinates. 
Also, we shall show how the iterative structure naturally associated 
with the (generalized) recursion relations defined by A^ allow to recur- 
sively construct the additional integrals of motion ensuring complete 
integrability. Finally, we will make contact with the "Lax" approach 
to the separability of the Neumann-Rosochatius system showing that 



(1.6) 



detiN{\a)) + PiXa) = 0, a = 1, . . . , n. 



(1.7) 



d{N*dH)(VH,VH) = 0. 
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the separation relations tying Darboux-Nijenhuis coordinates and these 
integrals are nothing but the spectral curve relations 

Obviously enough, the conditions on coming from ()1.7|) . in their 
full generality, are too difficult to be solved. The couple of Ansatze 
which will enable us to solve them for the NR case are the following. 
The first one is suggested from the fact that the phase space of the NR 
system is a cotangent bundle; accordingly, we will seek for a special ujN 
manifold structure on T*S^, defined by a (1, 1) tensor N induced by 
a suitable tensor (with zero torsion) L on the base manifold S"". The 
second one will be to use a special form of equation (|1.7p . that reads 

(1.8) d(N*dH - ^tr(Ar) A dH) = 0. 

The plan of the paper is as follows: In Section El we will collect some 
notions of the theory of ujN manifolds, and briefiy discuss the bi- 
Hamiltonian theorem for separation of variables. In Section |3] we will 
solve equations p.8|) . thus showing that the geometry of uN mani- 
folds can be used to discover the separation variables of the Neumann- 
Rosochatius (NR) system. In Section El we will find the Stackel sepa- 
ration relations and the family of commuting integrals of the system. 

Acknowledgments. We wish to thank Sergio Benenti, Franco Magri, 
and Giovanni Rastelli for useful discussions. The results presented in 
Section 2 have been obtained in collaboration with Franco Magri. This 
work has been partially supported by the Italian M.I.U.R. under the 
research project Geometry of Integrable Systems. 

2. The uN framework 

In this section we wish to recall some basic properties of a special 
class of bi-Hamiltonian manifolds, called utN manifolds. For a more 
detailed description we refer to |TBl lllj. By definition, an uN manifold 
is a smooth manifold M endowed with a pair of compatible Poisson 
bivectors Pq, Pi such that one of them (say, Pq) is nondegenerate. 
(Compatibility means that any linear combination of Pq and Pi is a 
Poisson bivector.) One can construct a recursion operator = PiPq^, 
whose Nijenhius torsion, 

(2.1) T{N){X, Y) = [NX, NY] - A^([A^X, Y] + [X, NY] - A^[X, Y]) , 

vanishes as a straightforward consequence of the compatibility between 
Pq and Pi (see, e.g., [14J). We set 2n = dim M and we denote by uq the 
symplectic structure associated to Pq, and by {■, - jo, {■, ■}i the Poisson 
brackets associated, respectively, to Pq, Pi. 
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The relevance of uoN manifolds in the theory of separable systems 
is mainly due to the existence, under suitable hypotheses, of a spe- 
cial class of canonical coordinates, that are selected by the geometric 
structure of the system itself. 

Definition 1. A system of local coordinates (xi, . . . , x„, yi, . . . , yn) that 
are canonical w.r.t. the symplectic form uq is said to be Darboux- 
Nijenhuis if the matrix expressing N in these coordinates is diagonal, 
i.e., 

iV = ^ {^i^ ® dxi + i^i-^ ® dyi^ . 

Notice that, since NPq is antisymmetric, it follows that Aj = z/j for all 
i. In general, however, the eigenvalues Ai, . . . , A„ need not be distinct. 

On the cotangent bundle of any differentiable manifold there is an 
elegant way to construct uN structures admitting Darboux-Nijenhuis 
coordinates through the following procedure [T2|. Let Q be an n- 
dimensional manifold equipped with a type (1, 1) tensor field L, whose 
Nijenhuis torsion vanishes. Let 6o be the Liouville 1-form and = d6o 
the standard symplectic 2-form on T*Q; the associated Poisson struc- 
ture will be denoted Pq. By thinking of L as an endomorphism of TQ, 
one can deform the Liouville 1-form to a 1-form 6^. 

{9l, Z)a = («, L(7r*Z))^(„) , 

for any vector field Z on T*Q and for any 1-form a on Q, where 
vr : T*Q ^ Q is the canonical projection. If we choose local coor- 
dinates (xi, . . . , Xn) on Q and set L{X) = Yl^j=i ^]-^^ 
local expression 6l = Yl^j=i L]yidxj w.r.t. the standard symplectic co- 
ordinates (xi, . . . , Xn, yi, . . . , yn) on T*Q. The complete lift of L is the 
endomorphism N of T{T*Q) uniquely determined by the condition 

(2.2) deL{X,Y)=uo{NXX), 

for all vector fields X, Y on T*Q. An easy computation shows that: 

(2.3) r,,A^.±n±-±y.i^-^,± 

II 

Since L has vanishing Nijenhuis torsion, the same property holds for 
the type (1,1) tensor field N on T*Q Prop. 5.6, p. 36], and 
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(T*Q, Pq, Pi := NPq) is an uN manifold [inj. The Poisson structure 
Pi is related to cji := dO^ by the formula: 

Pi((iF, dG) = uJiiXp, Xg) for all F, G G C°°(T*g) , 

where Xp, Xg are the Hamiltonian vector fields associated to F, G 
w.r.t. the symplectic form uq. By the very definition, if xj^"* is the 
Hamiltonian vector field associated to H w.r.t. Pi, then X^^ = NXh- 
Notice that, in general, this vector field need not be Hamiltonian or 
even locally Hamiltonian w.r.t. ujq. Indeed, the 1-form N*dH may fail 
to be closed (here A^* is the adjoint of the endomorphism A^), and one 
has 

(2.5) L^Xn^Q = -d{N*dH) = Lx„iui • 

In fact, from (Q it follows that {N*dF){Y) = -ujo{NXh,Y) = 
-ui{Xh,Y). 

Let us now assume that L has n functionally independent eigenvalues 
Ai,...,A„. Since L is torsionless, these eigenvalues determine local 
coordinates on Q satisfying the relations: 

dXi ^dXi 

We denote by /ij the conjugate momentum to A^; clearly, one has 
N-^ = Aj^. The coordinates (Ai, . . . , A„, fii, . . . , ^n) are Darboux- 
Nijenhuis coordinates for the uN manifold (T*Q, Pq, Pi). 

We can exploit the geometric setting of lvN manifolds in order to find 
an intrinsic separability condition for a given Hamiltonian function H G 
G°°{T*Q). Let us suppose that the vector fields Xh, NXh, ■ ■ ■ N'^'^Xh 
are pointwise linearly independent, so that they generate an n-dimensional 
distribution Vh. If we compute the conditions 

(2.6) d{N*dH){N'XH,N^XH) = for alH, j = 0, 1, . . . , n - 1 

in the Darboux-Nijenhuis coordinates (Ai, . . . , A„, fii, . . . , /i„), we get a 
system of differential equations equivalent to the Levi-Civita separa- 
bility formulae p. 208, eq. (1.230)]. 

Theorem 2. In the above hypotheses and notations, the Darboux- 
Nijenhuis coordinates associated with L are separation variables for H 
if and only if the 2-form d{N*dH) annihilates the distribution Vh- 

The separability condition ()2.6p implies that the distribution Vh is 
integrable. So, there exist n independent local functions Hi, . . . Hn that 
are constant on the leaves of Vh- The distribution being invariant un- 
der the action of N, the same is true for the differential ideal generated 
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by the HiS, so that the following condition holds: 



(2.7) 



N*dHi = y FikdH, 



k=l 



where is a matrix with distinct eigenvalues. Moreover, since Vh is 
Lagrangian with respect to both i^o we have: 



{Hi, H,}o = {Hi, H,}i = for all z, j 



It follows that the functions Hi, . . . H^ are separable in the the Darboux- 
Nijenhuis coordinates singled out by L. 

A particular case of this state of affairs is provided by the following 
example, which will turn out to be of great importance in the study of 
the NR system. Let us consider the characteristic polynomial 

(2.9) det(A/ - L) = A" - CiA""^ - CsA"-^ c„ 

of the endomorphism L, and assume we are given a Hamiltonian H 
satisfying the condition 

(2.10) d{N*dH) = dci AdH . 

From ()2.5j) it follows at once that this equation is equivalent to 

(2.11) d{LxjL- Hdci) = Q. 

The condition ()2.10|1 is a sufficient condition to the separability of H in 
the Darboux-Nijenhuis coordinates associated to N , because it implies 
that the 2-form d{N*dH) annihilates the distribution T>h generated by 
the vector fields Xh, NXh, ■ ■ ■ N'"'~^Xh- Moreover, it can be shown 
[THj that, choosing a local function H2 such that dH2 = N*dH — cidH , 
the 1-form N*dH2 — C2dH is again closed, so that we can find a local 
potential H^. By iterating this procedure, we end up with n indepen- 
dent local functions Hi = H, . . . , Hn that are constant on the leaves of 
T>H and satisfy the conditions: 



dHi+i = N*dHi - CidH 
(2.12) = N*dHr,- CndH . 

In this case the matrix F has the form 



C2 





1 





\Cn 



2,. 



0\ 



1 
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and the following condition is readily checked: 
(2.13) N*dF = FdF. 

We set F = diag(Ai, . . . , A„) 5*. Then, one has 5*-,^ = X^~^ and, 
by virtue of ll. Theorem 4.2], one obtains the separability equations 
'^k=i'^k^]~^ = Uj. Summing up, the functions Hi = H,...Hn are 
proved to be Stackel separable in the Darboux-Nijenhuis coordinates 
associated to L (this means that the separation relations are affine in 
the Hk's). 

3. TORSIONLESS TENSORS AND SEPARABILITY OF THE NR SYSTEM 

According to the results of Section 2, to separate the NR system 
we seek a tensor field L of type (1, 1) on 5"" satisfying the "strong" 
separability condition (j2.1Uj) . i.e., 

(3.1) d{LxjL-Hdci) = 0, 

where 9l = Ylab=i ^IPadqb (for any set of fibered coordinates {qa,Pa)) 
and Ci = trL. We have seen that the eigenvalues of such an L (if real 
and functionally independent) are separation variables for H. The form 
of the constraint and of the potential suggest using on S*" the coordi- 
nates Xa '■= Xa^, for a = 1, ... , n} If Ya are the momenta conjugated 
to the Xa (and the point particle has unit mass), the NR Hamiltonian 
is given hj H = T + V with 

T = 2Y,Xa{l-Xa)Yj -AY,XaXhYaY, 



V 



a<b 

(tta - an+l)Xa + 



/3n 



2 



+1 



2(1-Ea^a) 



Expanding in powers of the momenta, we see that condition 1)3.111 splits 
into 

(3.2) d{LxjL-Tdci) =0 

(3.3) d{Lxy9L - Vdci) = . 

Remark 3. As noticed in equation ()3.2)l means that L is a sym- 
metric conformal Killing tensor with respect to the usual Riemannian 
metric of 5", and implies that the torsion of L vanishes. On the other 
hand, equation ()3.3|1 can be written as d{L* dV + ci dV) = 0, which 
is a separability condition on the potential V appearing in the works 
of Benenti (see, e.g., [2J). However, since our approach applies also to 

^In this and in the foUowing section we use the following convention: middle 
indices like k run from 1 to ri + 1, while indices like a, b, c run from 1 to n. 
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systems which are defined on general symplectic manifolds (not nec- 
essarily cotangent bundles), or, in other words, to Hamiltonians that 
are not quadratic in the momenta, we will not use these results and 
we will solve directly equations ()3.2p and ()3.3p . We also observe that a 
significant part of the "Riemannian" theory of separation of variables 
can be seen as a particular case of the bi-Hamiltonian approach (see 
O E] and |SI , where the Neumann system is also discussed) . 

We start seeking a solution L whose dependence on the coordinates 
Xa is affine: = '^c^bc-^c + B^- Let us consider, for the sake of 
simplicity, the case n = 2. Condition ()3.2j) gives 

^12 ~ ^22 ~ ^11 ~ ^21 = = = , 
^21 — ^22 5 ^11 — ^12 — ^22 ^ B\ = B. 



2 

2 1 



SO that we are left with the unknowns A22, B\, and Now, condition 
()3.3p is equivalent to 

Al^{a2 - «i) = {Bl - B\){a^ - aa) , 

which means that = c{a^ — a2) and B2 — B\ = c{a2 — ai) for 
some constant c. Thus i?| = c(a2 + d) and Bj = c{ai + d), where d is 
another constant, and the components of L are given by 



t2 t2 
^2 



(as - ai)Xi + ai + d (0:3 - a2)Xi 

(as - ai)X2 (as - a2)X2 + a2 + d 



Since we are interested in the coordinates given by the eigenvalues of 
L, we can set c = 1 and ci = without loss of generality. 

Coming back to the general case, it is not difficult to check that the 
1-form 

)dXa , 

a \ b / a 

corresponding to the (1,1) tensor field given by 

(3.4) LI = (a„+i - ab)Xa + S^a^, 

satisfies conditions (j3.2p and (j3.3p . Although these formulas define L in 
coordinate patches, it is not difficult to show that L is globally defined 
on the whole S"'. Indeed, it is the restriction to of the tensor field 
L on M!^^^ defined as 

f d d Xi ^ 2 9 



dxi dxi ^ dxj 

j,k 
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where = Xi^. In order to show that L restricts to S*" it is sufficient 
to check that L*dr = 0, implying that, at every point of 5*", the image 
of L is (contained in) the tangent space to the sphere. 

Summarizing, we have found a tensor field L satisfying the separa- 
bility condition (jS.lj) : thanks to the result of ^ referred to in Remark 
01 the torsion of L vanishes. Thus the coordinates associated with L 
are separated variables for the NR system. Let us explicitly check that 
the eigenvalues of L coincide with the spheroconical coordinates. To 
this end we find convenient to introduce the following notations: let a 
and X denote the n-component vectors whose entries are, respectively, 

ctb = a„+i - ttfo, Xb = Xfo, 6 = 1,..., n, 

and let A be the n x n diagonal matrix of the parameters Then 
we can compactly write the matrix form ()3.4|) of the tensor field L as 

(3.5) L = A + X®ct. 

To compute the roots of det(A — L) we notice that the L is a rank 1 
perturbation of A; hence we write 

(3.6) A-L= (A-A) ■ (1-X'(A)®ck), 

where X'(A) is the vector with entries . Using the rank 1 

A — 

Aronszajn-Weinstein formula 

det(l + x®y) = 1 + (y,x), 

we arrive at 

(3.7) det(A -L) = n(A - «.)(! - i^^I^±lZ^). 

a=l 6=1 " 

Recalling the definitions Xf, = Xb^, for b = 1, . . . ,n, and the constraint 
Yli^i^ = 1; we can by means of elementary calculations conclude that 
such an equation is equivalent to 

(3.8) det(A -L) = - a,) T^' 
that is, the eigenvalues A^ of L satisfy the equations 

(3.9) Et^ = °' with^x^ = l. 



These are the well-known defining relations for the spheroconical (or 
elliptic-spherical) coordinates. 
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We close this section reporting, for the sake of completeness, the 
well-known computation of the momenta fia conjugated to the Aq. It 
is easily checked that the usual rule for computing the residues gives 



2 _ rial"* - ^a) 



(3.10) X, . 

Then we have that 

and, using again ()3.10|) . that 



1 \ ^ dXa 

2^' ^ Xn - 



If (xi, . . . , Xn+i,yu • • • , Vn+i) e M2n+2 3 TS^ ~ T*5", then the /x, are 
given by 

a \ i / \T*S" a \ i " V 



meaning that 



1 v;^ XiUi 



2 Aa - ttj 

Therefore, we can conclude that the separation variables Aa are the 
solutions of 

while the conjugated momenta are given by fia = h{Xa), with 

4. Integrals of motion and Stackel separability 

In the previous section we have found a tensor field L on 5" which 
gives the separation coordinates of the NR system (i.e., the spherocon- 
ical coordinates). Since L satisfies the "strong" separability condition 
(jH.lj) . we know from Section 2 that: 

(1) There is an iterative method for constructing n integrals of mo- 
tion in involution, {H = Hi, H2, . . . , Hn)- (Of course, we have 
to take into account that T*S"' is simply connected for n > 2.) 

(2) The NR system is Stackel-separable. 
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The integrals of motion are given by 
(4.1) dHa+l = N*dHa -CadH , a = 1, 



n 



where A" — Z]a=i CaA" " = det(AJ — L). This defines the Ha up to 
additive constants. For example, in the case n = 2 one finds 

H2 = 2{aiY2^X2 + a2Yi^Xi){Xi + X2 - 1) - 2^3X1X2(^2 - Yif 



a2(tti — «3jXi + ai{a2 — a3jX2 H — h 



Xi 



X2 



)X9 



!^1 



,02X1 + 0:1X2 



+ («3 - «2)/9l V + - «l)/92 V + 1 

Ai A2 i — Ai — A2 

Before showing that the Ha coincide with the integrals of motion known 
in the literature, let us consider the separability a la Stackel of the NR 
system. It is guaranteed from the results in Section 2 that 



(4.2) 



A^ ^Hb = Ua{\a, fJ'a) 



1 n 



6=1 



where the Ua are polynomials. 

Remark 4. It is easy to see that the polynomials X^ = Xa(Ai, . . . , A^j) 
and Ya = Ya{Xi, . . . , Xn, f^i, ■ ■ ■ , fJ^n) are invariant under the exchanges 
{Xb, fib) ^ (Xcfic)- This entails that Ub = Uc- 

Next we want to compare the constants of motion defined by ()4.1|1 
with the spectral invariants of the Lax matrix (see, e.g., [T^ ) 

-h{X) + ik{X) e(A) 
/(A) h{X)+ik{X) 



N{X) 



where h{X) = | has already been introduced, and 



X — a 



X — a,- 



The spectral invariants are the coefficients of the polynomial 

1 " 

(4.3) P(A) = a(A) det X(A) = tA" + ^ PaA""'' , 



a=l 



where a(A) = Hil-^ ~ ^i) restriction to T*S"' has been tacitly 

assumed. In particular. Pi = ^H, where H is the NR Hamiltonian. 
Our strategy to prove that Pa = \Ha for all a is to show that 

(4.4) N*dP{Xb) = XbdP{Xb) , 6 = 1, 



. n 
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which imphes that the Pa satisfy 

N*dPa = dPa+l +CadPi , a = 1, ... ,71 - 1 , 

because the A^'s are the roots of ()2.9p . Since these relations coincide 
with the equations ()4.H) for the Ha, and the starting points fulfill Pi = 
^Hi, we can conclude that 

(4.5) Pa = ^Ha ioT a = l,...,n. 

To show that (j4.4|) holds, we recall that fib = h{Xb) and f{Xb) = 0, so 

that ()4.3j) entails 

(4.6) 

P{Xb) = a(Afc) {-h{XbY - Khf - e{Xb)f{Xb)) = -a(A,) (/i,^ + k{Xbf) . 

Then (14. 4|) follows from the definition of DN coordinates. 

Finally, from ()4.5|) and ()4.6|) we obtain the separation relations for 
the Ha, 

n 

J2 HaXt" = -^K - MK) W + KXbf) , 

a=l 

i.e., the exphcit form of the Stackel vector with components Ub appear- 
ing in dOJ- 
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